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Abstract 

We study a singular limit for the compressible Navier-Stokes system when the Mach and Rossby 
numbers are proportional to certain powers of a small parameter e. If the Rossby number dominates 
the Mach number, the limit problem is represented by the 2-D incompressible Navier-Stokes system 
describing the horizontal motion of vertical averages of the velocity field. If they are of the same 
order then the limit problem turns out to be a linear, 2-D equation with a unique radially symmetric 
solution. The effect of the centrifugal force is taken into account. 

1 Introduction 

Rotating fluid systems appear in many applications of fluid mechanics, in particular in models of atmo- 
spheric and geophysical flows, see the monograph [3J. Earth's rotation, together with the influence of 
gravity and the fact that atmospheric Mach number is typically very small, give rise to a large variety 
of singular limit problems, where some of these characteristic numbers become large or tend to zero, see 
Klein [13]. We consider a simple situation, where the Rossby number is proportional to a small param- 
eter e, while the Mach number behaves like e m , with m > 1. Scaling of this type with various choices 
of m appears, for instance, in meteorological models (cf. [13j Section 1.3]). 

Neglecting the influence of the temperature we arrive at the following scaled Navier-Stokes system de- 
scribing the time evolution of the fluid density g = g(t, x) and the velocity field u = u(t, x): 

d t g + div x {gu) =0, (1.1) 

dt(gu) + dw x (gu ® u) + -(b x gu) + -^-V x p(g) = div a; S(Vxu) + -^gV x G, (1.2) 

e e e 
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where p is the pressure, and § is the viscous stress tensor determined by Newton's rheological law 

S(V x u) = /i(v x u + V^u - ^div^ul) , n > 0. (1.3) 

For the sake of simplicity, we have omitted possible influence of the so-called bulk viscosity component 
in the viscous stress. 

We consider a very simple geometry of the underlying physical space f2 C M 3 , namely f2 is an infinite 
slab, 

n = u 2 x (0,1). 

Moreover, to eliminate entirely the effect of the boundary on the motion, we prescribe the complete slip 
boundary conditions for the velocity field: 

u • n|en = 0, [Sn] x n\ an = 0, (1.4) 

where n denotes the outer normal vector to the boundary. Note that the more common no-slip boundary 
condition 

u|an = 

would yield a trivial result in the asymptotic limit, namely u — > for e — > 0. On the other hand, the 
so-called Navier's boundary condition 

u • n\ an = 0, ^u tan + [§n] tan | 9n = 0, /3 > 0, (1.5) 

gives rise to a friction term in the limit system known as Ekman's pumping, see Section [23] below. 

As is well known (see Ebin [7]), the boundary conditions (|1.4[) may be conveniently reformulated in terms 
of geometrical restrictions imposed on the state variables that arc periodic with respect to the vertical 
variable X3. More specifically, we take 

Q. = R 2 xT 1 , (1.6) 

where T 1 = [— 1, 1] | { 1. 1> i s a one-dimensional torus, on which the fluid density g as well as the horizontal 

component of the velocity = [u , u 2 ] are extended to be even in x 3 , while the vertical component u 3 
is taken odd: 

q(xi,x 2 ,-x 3 ) = g(xi,x 2 ,x 3 ), u l (x 1 ,x 2 ,-x 3 ) = u J (xi, x 2 , x 3 ), i = l,2, 

(1.7) 

u 3 (xi,x 2 , -x 3 ) = -u 3 (xi,x 2 ,x 3 ). 

Finally, we assume that the rotation axis is parallel to x 3 , namely b = [0, 0, 1], and set V X G ~ V^Ice^I 2 - 
the associated centrifugal force, where we have written Xh = [x 1 ,^ 2 ]. 

As shown in [3] , incompressible rotating fluids stabilize to a 2D motion described by the vertical averages 
of the velocity provided the Rossby number e is small enough. Besides, the stabilizing effect of rotation 
has been exploited by many authors, see e.g. Babin, Mahalov and Nicolacnko [T], [2]. On the other 
hand, compressible fluid flows in the low Mach number regime behave like the incompressible ones, see 
Klainerman and Majda [12] , Lions and Masmoudi [17] , among many others. Thus, at at least for m 1, 
solutions of the scaled system (|l.ll) . (|1.2[) are first rapidly driven to incompressibility and then stabilize 
to a purely horizontal motion as e — > 0. On the other hand, the above mentioned scenario changes 
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completely if m = 1. In this case the speed of rotation and incompressibility act on the same scale. 
Accordingly, the limit behaviour of the fluid is described by a single (linear) equation, see Section [2~5l 
Note that such a picture is in sharp contrast with [3] , where the effect of the centrifugal force is neglected. 

However, a rigorous justification of the above programme is hampered by serious mathematical difficulties: 

• The main issue in the low Mach number limit, at least in the case of the so-called ill prepared initial 
data, is the presence of rapidly oscillating acoustic waves, cf. Desjardins and Grenier [5], Desjardins 
et al. [6], Lions and Masmoudi [17]. Similarly to [5], given the geometry of the spatial domain f2, 
we may expect a local decay of the acoustic energy as a result of dispersive effects. Unfortunately, 
however, the fluid is driven by the centrifugal force that becomes large for \xh\ — > oo. Specifically, 
we have G « e~ 2m on the sphere of the radius e~ m , whereas the speed of sound in the fluid 
is proportional to e~ m . In other words, the centrifugal force changes effectively propagation of 
acoustic waves and this effect cannot be neglected, not even on compact subsets of the physical 
domain. 

• The dispersive estimates of Strichartz' type exploited in [3] cannot be used in the present setting 
as the acoustic waves, represented by the gradient component of the velocity field, remain large 
for \xh\ — > oo. 

Our approach is based on combination of dispersive estimates for acoustic waves with the local method 
developed in |llj . As already mentioned, we focus on two qualitatively different situations: one where 
m>l, and one where m = 1. 

• The case m > 1. 

In order to eliminate the effect of the centrifugal force, we compute the exact rate of the local decay 
of acoustic energy, here proportional to e m / 2 , by adapting the argument of Metcalfe [20] (cf. also 
D'Ancona and Racke [4], Smith and Sogge [21] ). Accordingly, the associated acoustic equation can 
be localized to balls of radii e~ a for a certain a > m/2, on which G(x) w e~ 2a . Having established 
local decay of acoustic waves, we use the method developed in [TT], based on cancellation of several 
quantities in the convective term, similar to the local approach of Lions and Masmoudi [18] . 

• The case m = 1. 

In this case both the high rotation and weak compressibility limits occur at the same scale. We can 
compute immediately the limiting diagnostic equations, and, similarly to the previous case, careful 
analysis of cancellations in the convective term allows to identifiy the limit system. In contrast 
with the situation studied in [9] , the limit system is linear as a consequence of strong stratification 
caused by the centrifugal force. 

The paper is organized in the following way. In Section [21 we introduce the concept of finite energy 
weak solutions to the compressible Navier-Stokes system (jl.ll - [L4)l and recall their basic properties. In 
particular, we present uniform bounds on solutions independent of the scaling parameter e — > 0. The main 
results concerning the singular limit of solutions to (jl.ll - ll.4j) are stated in Section [2] The anisotropic 
situation, when m 3> 1, is analyzed in Section [3] by means of several steps: The easy part concerns 
the formal identification of the limit system, while, in Section 13.21 the propagation of acoustic waves is 
studied as well as their local decay in f2. Finally, the limit systems for m 3> 1 is justified in Section [3. 31 by 
means of a careful analysis of the convective term. Finally, the isotropic case, when m = 1, is examined 
in Section [4[ where the proofs borrow several ingredients introduced in jllj . 
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2 Preliminaries and statement of the main results 



In this section we introduce the main hypotheses, recall some known results concerning existence of 
solutions to the primitive system as well as the nowadays standard uniform bounds independent of the 
scaling parameter, and, finally, formulate our main results. 

2.1 Main assumptions 

Consider a family of solutions g e , u e of the Navier-Stokes system - II .4[) in (0, T) x 0, emanating 
from the initial data 

Qe(0,-) = Q0,e, U E (0,-) = U 0)£ . 

We assume that the initial data arc ill-prepared, specifically, 

eo,e = Q e +e m r . e , (2.1) 
where g £ is a solution of the associated static problem: 

V x p(&) = £ 2(m_1) &V x G in 0. 

Consequently, 

P(q c ) = e 2{m - 1} G + const, where P(g) = f dz. 

Ji * 

Furthermore, we suppose that 

p G C^OjOo) nG 2 (0,oo), p'(g) > for all g> 0, lim = c> (2.2) 

for a certain 7 > 1 specified below, and that 

Ge^(n), G(x)>0, G(x 1 ,x 2 ,-x 3 ) = G(x 1 ,x 2 ,x 3 ), 
|V x G(.t)| < c(l + |ar fc |) for all x G 0. 

Finally, we normalize 



(2.3) 



P(£ e ) = e 2 (™- 1 )G, (2.4) 
noticing that 

Qe{x) > 1, Qe{x) — > 1 for any x & Q as e — > for m > 1, (2.5) 

whereas 

f} £ = p is independent of e provided m = 1. (2-6) 
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2.2 Energy inequality 

Introducing the relative entropy 

E(g, g) := H(g) - H'(g)(g - g) - H(g), H(g) := g f ^-dz, 

Ji z 

we note that H'(g) = P(g) + const; therefore we may assume that 

jf Qfo|u e | 2 + -^E( Qe , g £ )j (t,-) + ^ ^S(V*u £ ) : V x u £ dxdr' 

< Q<?o, £ |u , £ | 2 + -^E(q q>£ , £ £ )^ dx 
including implicitly the mass compatibility condition 

(Qs ~ Qe) d.T = 0. 



(2.7) 



Now, in order to establish uniform bounds independent of the scaling parameter e, the initial data 
specified in (|2.ip must be chosen in such a way that the expression on the right-hand side of (|2.7[) remains 
bounded uniformly for e — > 0. Thus, if 7 < 2 in (|2.2[) , it is enough to assume that 

{ro, E }e>o bounded in I? n {\/^e u o,e}£>o bounded in L 2 (f2; i? 3 ). 

In general, we suppose that 

\g~ e ~ L 2~ro,\ bounded in L 2 (Vl), {r „} £ >o bounded in L 2 n L°°(Q), 

l ' J £>0 

{v / ge"uo, e } £ >o bounded in L 2 (il;R 3 ). 



(2.8) 



2.3 Finite energy weak solutions 

We say that g, u is a finite energy weak solution of the Navier-Stokes system (11.11 - [T31) in (0, T) x SI, 
supplemented with the initial data (|2.1[) . if: 

• the energy inequality 

f a (^l u l 2 + ^L E (e> &)) ( T > ') dx + f J n § ( v * u ) : v * u dx dr ' 

< ^£>o, £ |u , £ | 2 + ■^E(g 0>e , g e )\ dx 

holds for a.a. r e (0,T); 

• equation (|1.1[) is satisfied in the sense of distributions, specifically, 

/ / (gd t (p + gu- V x (p) dx dt = - / go, e ^(0, •) dx 
for any V3 <E C* c °°([0,T) x 0); 
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• the pressure p is locally integrable in [0, T) x 0, equation (|1.2j) holds in the sense of distributions: 
/ / (gu ■ d t ip + gu (g) u : V x </3 — -b x (gu) • ip H — ^— p(g)div x ip) dx dt 

Jo y x £ 7 (2-9) 

= / / (§(V x u) : V x <p -gS/ x G-ip) dx dt - / £>o,eU , e • tp(0, •) da; 

Jo Jn K £ ' Jn 

for any ip € C* c °°([0,T) x 0; i? 3 ). 

Under hypothesis (|2.2p . the existence of finite energy weak solutions can be established by the method 
developed by Lions [Hj , with the necessary modifications specified in [Til] in order to handle the physically 
relevant range of adiabatic exponents 7 > 3/2. 



2.4 Uniform bounds 

In order to study the asymptotic behavior of solutions, we first establish uniform bounds independent of 
the scaling parameter £ — > 0. As a matter of fact, all of them follow from the energy inequality ()2.7|) . 
Similarly to Chapter 5] , we introduce the essential and residual component of a function h as 



h = h P 



h 

"'res: 



where 



h ess (t,x) = h(t,x) for (t,x) such that g E (t,x) £ (1/2,2), h css (t,x) = otherwise. 
Now, by virtue of (|2~4"1) . 

1 < Qe{x) < 1 + c(r)£ 2(m ~ 1 ~ Q) for all x e B r/£a , < a < m - 1, 

where we have denoted 

B R = {xen I |o; A | <R}. 
It follows directly from energy inequality (|2.7[) that when m > 1 + a 



ess sup 

te(o,T) 





& - Qe 








CSS 



ess sup 

te(0,T)JB r/E 



ess sup 
te(o,T) 



< c(r) 
[fe]? es dx < £ 2m c(r) 
l ros dx < e 2m c(r) 



and 



ess sup |{x e I e (i,a;) < 1/2}| < ce 2m < c 
te(o,T) 

ess Sup llv^Uelli 2 ^;^ 3 ) < c 

te(o,T) 



/ / §(V x u £ ) : V x u £ dx dt < c. 
Jo Jn 



(2.10) 



(2.11) 
(2.12) 

(2.13) 
(2.14) 
(2.15) 

(2.16) 



(> 



In the case when m = 1, the bounds (|2.1ip - (|2.13p should be replaced by 

< c(r) 



ess sup 
te(o,T) 





' Qe - Pe" 








ess 



CSS 



sup f [g e }l s dx < e 2m c(r) 

i(0,T) JB r 

SUp / l res dx < £ 2 ™c(r). 
E(0,T) ifl r 



te(0,T) JB r 

ess 



(2.17) 
(2.18) 
(2.19) 



te(o,T) Js r 

Finally, combining (|2.14l - l2~.16p with a variant of Korn's inequality (see [HI Theorem 10.17]), we obtain 

r-T 



U 



di < 



(2.20) 



All generic constants in the previous estimates are independent of the scaling parameter e. 
2.5 Main results 

In what follows, the symbol H denotes the standard Hclmholtz projection onto the space of solcnoidal 
functions in O, specifically, 



H[v]fcfc) = v(£fc) 



1 



e(* • fc) + fc)) , *(e ■ v h (£, fc) + fc« 3 (e, fc) 



H^fvl = v - Hfvl 



where the symbol v(£,fc), £ S i? 2 ,fc £ Z, denotes the Fourier transform of v = v(xh,X3). Similarly, the 
Laplace operator A is identified through 

Kv « -(|e| 2 + fc 2 ) v(£,k). 

Finally, we introduce the vertical average of a function v as 

1 



(v) (xh) 



IT 1 ! 



v(x h ,x 3 ) dx 3 . 



(2.21) 



T 



2.5.1 Multiscale limit (m > 1) 



Theorem 1 Let the pressure p and the potential of the driving force G satisfy hypotheses h2.2\) . i2.3\) . 
with 7 > 3/2. Let g £ , u £ be a finite energy weak solution of the Navier-Stokes system in (0,T) x 
belonging to the symmetry class jl.7\) , emanating from the initial data h2.1\) . 12.8]) . In addition, suppose 
that 

m > 10 



and that 
Then 



Uo j£ —> Uq weakly in L 2 (fl;R 3 ). 



ess sup \\g e — lj| (L 2 +L-t)(K) < £ m c{K) for any compact K C 0, 
te(o,T) 
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u e -> U wea/c/y in L 2 (0, T; W 1,2 (fi; i? 3 )), 
where U = [Uft(xh), 0] is </ie unique solution to the 2D incompressible Navier- Stokes system 

div^U,, = 0, (2.22) 

dtUh + dw h {U h ® U h ) + V h n = /iAftU/,, (2.23) 

loitft i/ie initial data 

U h (0,-) = [H[[(U ) h) 0] 

Remark l2.5H ^4 s/iort inspection of the proof of Theorem^ given in Section^ below reveals that replacing 
the complete slip condition \1.J$ by the Navier 's slip condition M.5\) would produce an extra term /3TJh 
on the left-hand side of (2.23\) known as Ekman's pumping. 

2.5.2 Stratified limit (to = 1) 

Theorem 2 Let the pressure p satisfy hypotheses (2.2\) , with 7 > 3 and let G(xf l ) = |x^| 2 . Let g s , u e 
be a finite energy weak solution of the Navier-Stokes system in (0, T) X Q belonging to the symmetry 
class \1. 7p , emanating from the initial data i2.1\) , h2.8\) , where 

to = 1. 

Ln addition, suppose that 

ro, e — > ro weakly in L 2 (f2), Uo >e — > Uo weakly in L 2 (Sl;R 3 ). 

Then 



r weakly-(*) in L°°(0, T; L 2 (K)) for any compact K C f2, 

t 

u e -> U weafcfy in L 2 (0, T; W 1,2 (fi; i? 3 )), 

wii/i 

?* = r(t,Xh) radially symmetric , U = [\Jh(t, Xh), 0], (2.24) 

and 

V ft (P'(e)r)+U^ = 0. (2.25) 

Moreover, the function r satisfies 

d t (r-div h (gV h (P'(Q)r))) +A 2 (p'(g)r) =0. (2.26) 
In addition, the initial value r(0) is i/ie unique radially symmetric function satisfying the integral identity 

j (gVfc (P'(e)r(0)) ■ + r(0) V>) d.T h . = jf^ ( (gUo.fc) ■ V^V + (r ) V>) dx ft (2.27) 

/or radially symmetric ip — ij)(xh) £ C^°(i? 2 ). 
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Remark 12.51 2 The initial condition (|2.27[) can be interpreted in polar coordinates as 




where s — \xh\, provided all quantities are sufficiently smooth. 

Remark [2TB1 3 Note that V2J>ty impl ies that V ' ug ■ = 0, or, equivalently, 

U/j • Xh = 0, meaning the limit velocity U/,. is tangent to the level sets {\xh\ = const}, 

Moreover it follows from \ 6 2,.2J$}ll3^ that |U/j, | is constant on \\xh\ = const}, and div^U^ = 0. 

It is remarkable that the limit equation (|2.26p is linear for m = 1, in sharp contrast with the homogeneous 
case g = const treated in [9]. This is related to the fact that the limit density g is stratified (non- 
constant). More precisely, the absence of nonlinearity is related to the smallness of the kernel of the 
penalized operator, defined by (|2.25[) . This phenomenon had already been identified for fluids with 
variable rotation axis, see [TT]. Much of the analysis that we shall follow to prove Theorem [3] borrows to 
this last reference. 

The rest of the paper is devoted to the proof of Theorems Q] and [U in Sections [3] and @] respectively 

3 Anisotropic scaling: Proof of Theorem [I] 
3.1 Preliminary remarks 

We start with some simple observations that follow directly from the uniform bounds (|2.11l - I2.20[) . 
Clearly, relations (|2~T0|) , (|2~TTj) . and (|2~T2)) imply that 

ess sup \\g E — 1\\(l 2 +li)(k) < e m c(K) for any compact K C 0, 
te(o,T) 

while (|2.20p yields immediately 

u e — > U weakly in £ 2 (0, T; W 1,2 (fl; R 3 )), at least for a suitable subsequence. 

Consequently, letting e — > in yields 

div x U = a.a. in (0, T) x O. (3.1) 

Moreover, it follows from equation (|1.2p that 

H[b x U] = 0, meaning, b x U = V x $ for a certain potential $. 

Consequently, $ and arc independent of the vertical coordinate £3, and moreover, divhXJh = 0. 
Since U is solenoidal, we have d X3 U 3 = 0; whence, as U 3 has zero vertical mean, 

U 3 = 0, U=[V h (t,x h ),0]. 

Finally, in view of the uniform bounds established in (|2.11[) , (|2.20p . it is easy to pass to the limit in 
all terms appearing in the weak formulation of momentum equation (|2.9p . tested on ip = [4>(t, Xh), 0], 
with div/j</> = 0, with the exception of the convective term div x (g e u e u £ ) that will be analyzed in the 
next two sections. 
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1 + ~ <<*< \{m-2)- (3.2) 



3.2 Acoustic waves 

Since m > 10, there exists a such that 



This choice of the parameter a will become clear in Section 13.31 Moreover, we introduce a family of 
cut-off functions Xe such that 

Xe G CT(i? 2 ), < Xe < 1, Xe(z h ) = 1 for |^| < -1 , 

£ (3.3) 

X e (x h ) = for \x h \ > A |V x x e (a:fc)| < 2e Q for s fc G i? 2 . 

As the density becomes constant in the asymptotic limit, the basic idea is to "replace" u £ « x £ p £ u £ and 
write 

XeQe^e = H[x E £ £ U £ ] + ^[XeQeUe] = (H[x £ £> £ U £ ]) + (H[x £ g £ U £ ] - (H[x £ £> £ U £ ]) ) + V s # £ , 

where V K * e = H J -[x £ £' e u e ]. 

In the remaining part of this section, we examine the asymptotic behavior of the acoustic potential ty e . 
More specifically, we show that Vi;$ E tends to zero on compact subsets of f2 and therefore becomes 
negligible in the asymptotic limit s 0. 

3.2.1 Acoustic equation 

Following Lighthill [II], [15], we rewrite the Navier-Stokes system (jl.ip . (|1.2[) in the form: 

£ m^&_l£^ + div,(£, £ U £ )=0, 

£ m 5 t ( & u £ )+/(i)v x (^^) =- i ^v ;c (p( & )-p(fc)-P , (i)(es-^)) 

+ e m di V:c (§(V K u £ ) - ( fc u £ x u £ )) - e" 1 - 1 (b x & u E ) + e^ m ^^^V x G 

where, exactly as in Section 12.31 equations are understood in the sense of distributions in (0, T) x O. 
Furthermore, introducing new variables 

a Qe - Qe 

be = Xe — , m e = XeQeU E , 

where Xe is the cut-off function specified in (|3.3p , we arrive at the equation 

e m d t S £ + div x m e = V xX e ■ (&Ue), (3.4) 

while m £ satisfies 

e m d t m £ +p'(l)V x S s = -e™- 1 (b x m e ) + ^"^^V.G + P '(l)V xX J e ~ & 



e m div x {xeS(V x U e ) - Xe(QeU £ X U £ )) - £ m (§{V x U e ) - Q s U e X U £ ) • V xX e (3.5) 
Xe(p(Qe) -p(Qe) -Jj'(l)(& ~ Qe)) + —{l>(.Qe) ~ P(t>e) ~ Jj'( 1 )(& - 0e))V a! Xe- 
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3.2.2 Uniform bounds 

Our next goal is to deduce uniform bounds on all quantities appearing in the acoustic equation 
([33]) . To begin, it follows from (|27TTT - 1235]) that 

{S £ } £>0 is bounded in L°°(0, T; L 2 + L 1 (0)). 

As a matter of fact, we have 

{5 e } e>0 is bounded in L°°(0,T; (L 2 + L 1 n L 7 )(^)) 
as the "residual set" is of small measure, cf. (12.131) . 
Similarly, we deduce from (j2~T2"j) . ([2~T5]) that 

27 



(3.6) 



{m e } E>0 is bounded in i°°(0, T; i + L g (D,; R )), q 
Moreover, combining (|3.3[) . hypothesis (|2 .3[) and (|3 .6[) we may infer that 



7 + 1 



V^G^ is bounded in L°°(0,T;L Z + L 1 ^)) 



;>0 



Finally, by virtue of (|2~T5l) . (|2~^0j) . 

{S(V x u £ )} E>0 is bounded in L 2 (0, T; L 2 (fi; i? 3x3 )), 

and 



{g E u e <8 u e } e >o is bounded in i°°(0, T; ^(Q; i? 3x3 )). 
Now let us estimate the pressure perturbation. Writing 

P(Qs) ~P(Qs) - ffe) =P(&) -p'(Qe)(Qe - Qs) + (p'{Qs) -p'(l))(Qs - Qs) 

we have 

P(Qs)-p(Qe)-p'{ee)(Qe- Qe) = p{Qs) ~ p(Qs) ~ p' '(Qs){Qs ~ Qs) + p(Qs) - p(§s) - P - Qs) 

J ess 

where, since p is twice continuously differentiable in (0,oo) and g e satisfies (|2.10p . 

p(Qs) ~ p{Qe) ~ p'{Qe){Qe ~ Qs) (t, x) < c [g e - g £ ] 2 css (t, x) provided x £ B 2/sa . 

L J ess 

Similarly, by virtue of (|2~T2"j) . (|2~T5)) 



(3.7) 

(3.8) 

(3.9) 
(3.10) 



ess sup 

te(o,T) Jb 2/e 



P(Qs) -P(Qs) -p'i.Qe){Qe - Qs) 



dx < ce 



Finally, in accordance with (|2.10l) , 

\p'(Qs(x)) -p'(l)\ < c\q e (x) - 1| < ce 2 ^- 1 -") for x e B 2/e a 
Thus, summing up the previous estimates, we conclude that 



ess sup 

te(o,T) 



^(p(Qs) -p(Qs) -p'{1)(q s - Qe)) 



<c(l + e m ~ 2 ~ 2a ). (3.11) 



(Li+L2)(B 2/eQ ;fl3) 
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3.2.3 Spatial regularization 

In view of the uniform bounds obtained in the previous section, equations (|3.4[) . Q3.5P can be written in 
the form 

e m d t S E + dW x m e =e a F}, (3.12) 
e m d t m E + P '{1)V X S E = (e m + e^-^JdiVxFjj + (e m_1 + e a + £ 2 ( m_1 ~ Q ))F^, (3.13) 

with 

{F*} e>0 bounded in L°°(0,T; L 2 + L x (0)) 
{F 2 } e>0 bounded in L 2 (0, T; L 2 + L 1 (f2; i? 3x3 )) 

{F 3 } £>0 bounded in L 2 (0, T; L 2 + L 1 (fi; R 3 ))- 

Equations p,12j) . p,13p are satisfied in the sense of distributions. For future analysis, however, it is more 
convenient to deal with classical (smooth) solutions. To this end, we introduce 

vs = K.S * v, 

where Kg = k$(x) is a family of regularizing kernels acting in the spatial variable. Accordingly, we may 
regularize (|3. 12|) . (|3.13[) to obtain 



s m d t S £i5 + div x m £ , s = e a Fl s 
e m d t m eiS +p'(l)V x S E ,s = (e m + e 2(m - 1 - Q) )div ;E F^ + (e m_1 + e a + £ 2 ( Tn_1_Q ))F 3 s 



where 



F l8\\LH0,T;W<°-*(a) < c(k,6) 



\Wl.s\\L 2 (0,T-W k . 2 (a:R 3x3 )) < c (k,S), 
\\ F3 eJ\L 2 (0,T-W k - 2 (n;R 3 ) < c(k,S) j 

for any k = 0, 1, . . . uniformly for e — > 0. 



(3.14) 
(3.15) 

(3.16) 



3.2.4 New formulation of the regularized acoustic wave equation 

Our aim is to rewrite the regularized acoustic equation (|3.14[) . (|3.15[) in terms of the acoustic potential. 
To this end, we decompose 

m £ ,«5 =Y e ,S + VxV E ,6, 



where Y 



e,S 



Hm 



e,<5j 



Introducing new unknowns S St s, ^e.s, we rnay rewrite the acoustic equation (|3.14[) . (|3.15[) in the form: 

e m d t S E ,s + A^ e j=e a Fl 5 , (3.17) 

e m d t * eiS +p'(l)S e ,s = (£ m +e 2 ^- 1 - a ))A-Miv x div x [F^]+( e ro - 1 + £ a + e 2 (" 1 - 1 - a ))A- 1 div x F^. (3.18) 

Our aim is to use dispersive estimates for (|3.1T[) . (|3. 18|) to deduce local decay of the acoustic potential. 
To this end, we present a simple result, the proof of which is a straightforward adaptation of Metcalfe [20l 
Lemma 4.1] (cf. also D'Ancona and Racke [U Example 1.2]) : 
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Lemma 3.1 Consider <p £ C c °°(i? 2 ). Then 

<p(x h ) exp (iy/^Atj [v] dx dt < c(^)||u||| 2(0) . 



oo JQ 



(3.19) 



Proof: For a function w = w(t,Xh,X3), wc denote by u>(r, £, k) its Fourier transform in the space-time 
variables, t £ R, !; £ R 2 , k £ Z. Accordingly, by virtue of Parseval's identity, 



oo J Q 



ip{x h ) exp (iV-At) [v] 



dx dt 



; E 
: E 



R 2 



— ? 7)^( T — Vlvl 2 + k 2 )v(r], k) dr] 



d£ dr 



oo -/i? 2 



{r=\/hl 2 +fe 2 } 



d£ dr. 



Furthermore, by the Cauchy-Schwartz inequality, 



E 



kez J -°° 
< 



R' 2 



- ri)v{r),k) dSn 



dC dr 



1^-77)1 dS v 



E 

ez 

we/ r i , 

(¥>)£ / / |^- ? /)||^,fc)| 2 dr ? de< C M||«|li 2( o)- 



^E L /°° / „_„,,„ i^-^)ii^ fc )i 2 

tez' 
fcez ' 

That proves Lemma 13. II 



m-v)\\v(r,,k)\ 2 dSjdZdT 



Q.E.D. 



The gradient of the acoustic potential <F e ,,5 can be written by means of Duhamcl's formula: 



1 ( ■ 1 
-expl-iV-A— 



-t - s 



exp iV-A— — - exp -iV-A 



-t - s 



ds 



1 + e 



m-1-1a rt 



-t - S 



exp iv— A- 

V \ t 



exp — iv— A- 



-t-s 



-1 _j_ g.CK— 7T1 _|_ p77l— 2— 2a pt 



-P^V-A^j+exp^-iV-A £ 



[V x A- 1 div a .div a .I^ ) j] ds 
t- s s 



[V^A-Miv^F^ s ] ds 
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where, for the sake of simplicity, we have set p'(l) = 1. Now, in accordance with Lemma \'3.1l 



10 JK 

and, similarly, 



exp ( iV-A— ) [ v ] 



dx dt < e r 







Ik 


L exp (* 










< T 1 




J 


3 JO JK 



JK 



-t - S 



exp (h/^At ) [v] 



dx dt < e r ' 



c v 7 2 



L 2 (fi)' 



(3.20) 



exp I A- 



[ 5 ( S )] ds 
-t - s 



dx dt 



[<?(«)] 



da; dt ds 



(3.21) 



<c7V 



exp 



(-if) 



L 2 (0) 



£m ||ffllL 2 ((0,T)xf2) 



for any compact A" C fi. Combining (|3.20[) . (|3 . 2 1 1) with the uniform bounds (|3.16[) and hypothe- 
ses flUJ, dUBl we infer that 

r T 

/ \\V x y £ ,6\\h(K-R3) dt < e 2/3 c(<5, A, T) for any compact A" C fi (3.22) 
Jo 

uniformly for e — > 0, where /3 > 1 provided a satisfies (|3.2[) . Thus the effect of acoustic waves becomes 
negligible in the limit e — > 0. Accordingly, the information about the asymptotic behavior is provided 
by the solenoidal component of the velocity field analyzed in the following section. 

3.3 Solenoidal part 

In order to control the solenoidal component of the velocity field, we write the momentum equation (|3.5p 
in the form: 



ed t m e $ + b x m e g = sdiv x 



Xe — V,G 



Xe§(V 2; U e ) - XeiQe^e X U e ) - £ (S(V X U £ ) - Q s ll e X U £ J ■ V^Xe 
l ^ X [Xe(p(Qe)-p(ge))}s - [V xX s(p(Qs) ~ p(Qs))}5 



in other words 

e&m e ,* + bxin ei tf=(e + e m - 1 - a )Qe, i +e 1 ^ (3.23) 
where 

{Qe,4e>o is bounded in L 2 (0, T; W k ' 2 (Q; A 3 )) for any fixed k, 5 > 0. 

3.3.1 Compactness of vertical averages 

Taking the vertical average of equation f|3 . 23[) we obtain 

ed t (m M ) + b x (m £ , 5 ) = (e + s™- 1 -") (Q e><5 ) 

+ e 1 - 2m fv :c ([x £ (p(&) -P(e=))]«) - ([Vxx £ (p(&) -p(e e ))]«> 



(3.24) 
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Recalling 

we check easily that 



m e S = Y s>s + V x # e ,a, Y e><5 = H[m £ , 5 ] 



bx (Y 



-(Ye* 8 



Y 1 

X £,8 







is an exact gradient as Y £j( 5 is solenoidal. 

Consequently, testing equation (|3.24p on ip 6 C£°(Q; i? 3 ), div x ip = we get 

d t [ (m e , 5 ) ■ tp dx = I (Q e ,j) • dx - - / (bx V. t f E i) ■ ^ da; 

provided e is small enough so that Xelsupp^ = 1- Thus we may use (|3.22l) to conclude that 

(Y Ej s) U5 strongly in L 2 ((0, T) x K; R 2 ) for any compact K C il and any fixed S > 0. (3.25) 
We note that this step depends essentially on (|3.2j) that requires the rather strong assumption m > 10. 



3.3.2 Oscillations 

We write any function v in the form 

v(x) = (v) (x h ) + {v} (x). 
Since {v} = v — (v) has zero vertical mean, it can be written in the form 



{v} (x) = d X3 I[v] with / I[v](x) dx 3 = 0. 

Jr 1 



(3.26) 



Moreover, we define 



u e,5 



Going back to equation (|3.23|) . we deduce that 



edtvlj + div h [Y EjS ]h = e[d Xl Ql s - d X2 Ql s ) - A h V E , 5 



and, finally, 
for all t e (0,T), x g B 1/e 



edt^J + d X3 Y 2 s = e{d Xl Q% s - d X3 Ql s ) - % 3tWa % 
ed t u 2 J - d X3 Y e )s = e(d X2 Ql 5 - i), ,(/-'.,) + % aiXl * 



(3.27) 
(3.28) 

(3.29) 
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3.3.3 Analysis of the convective term 

Following step by step the analysis performed in [TTJ Section 3] we observe that the only problematic 
component of the convective term reads 



H 



• Step 1: 

Since Y e $ is solenoidal, we have 



/ [div x (Y E>s ^Y e ,s)] h dx 3 
Jr 1 



div x (Y e , 5 <g> Y e>s ) = -V^Y^I 2 - Y £i5 x (curl[Y £)5 ]). 
As the former term is a gradient, we concentrate on the latter. 

• Step 2: 

Write 

Y £t5 x curl[Y E , 4 ] = (Y e ,j) x curl (Y e , 5 ) + d X3 ( (Y e , s ) x curl I[Y e , 5 ] + I\Y S>S ] x curl (Y s 
+ d X3 I[Y £>s ] x c^curl^Y^]), 



where the term in the brackets has zero vertical mean, while, in accordance with (|3.25|) . the first term is 
pre-compact. Finally, we have 



d X3 I[Y £)S ] x d X3 cvLvl(I[Y e , 5 }) 
= d X3 i[Yl s ]d X3 (d Xi i[Yi s } - i), / ) ".„ ) = d X3 I[Yl s }d X3 I[J e % j = 1,2,3. 



• Step 3: 

In accordance with p. 271 - 13.29[) , we get 



ed t {d X3 I[u>l>*}) + dl 3 I[Yl 5 ] = e(d Xl {Ql. s - (Qh)) - d X3 Ql, s ) - 



and 



ed t (d X3 I[u 2 J}) - d 2 X3 I[Y e ] s ] = e(d X2 (Ql s - (Q* 4 » - d X3 Ql s ) - ^„ ai * e>4 
at least for x e Bj^n. Next, compute 

[d X3 I[Y etS ] x c^curl^Y^])] 1 = d x J[Yl s ]d X3 I[u:lj] + d X3 I[Y E 3 s ]d X3 I[u 3 e '}] 

r2 1 rr, ,2,11^ p.2 T\V 2 lT\,, 2 'h rrj;.. rv .1 la rr. .3,li 



d X3 [d X3 I[Y* s )I[^ s ]) - d 2 X3 I[Y^ 5 }I{u^} - I[div h \Y e , s ] h ]d X3 I[u,%], 
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where, furthermore, 

^,J[i^,]/[<' 4 1 ] + /[div h [y e ,,] /i ]a a 3j[ w y] = - £ s t (9 l3 /[^ 3 ])/[ WE 2 j 



3,li 

s\ 



-si 



2,1-. 



for t G (0,T), x G B\j e a. Treating the second component in a similar fashion, we conclude that 

ttptt / / [div :c (Y E)5 <g> Y e ,j)] ■ <p dx dt -> / / div^Ua <£> Us) • dx h dt as e -> 
I' Ijo Jn Jo Ji? 2 

for any fixed 5 > 0, and for any 93 = £/>), 0], G C~([0, T) x A 2 ; A 2 ), div/^ = 0. 



(3.30) 



3.4 Convergence - conclusion 

In order to complete the proof of Theorem [TJ write 

Q e U e (g) U e = (£> £ - l)u e (g) U £ + (u £ j - Ug) (8) U e> 5 + U e ,<5 <g> (u e>( 5 - U e ) 

+ [(1 - e e )us]5 <X> u e>4 + u Ej5 (8 [(1 - g e )\i e ]sm £! s ® m e S 
for t G (0, T) and x G A' C O, A compact. 
Now, it is easy to observe that 

(g e - l)u £ <g u e , [(1 - £» e )u e ] 5 ® u £i5 -> in L 1 ((0, T) x A") for e -> 

while 

||u £i(5 - UsIlia^jBS) < cd\\u e \\ w i,2(K;R?) uniformly for e > 0. 
Consequently, testing the momentum equation ()1.2|) on a compactly supported solenoidal function, the 
convective term g E u E (g> u e can be replaced by m e ,5 <g> m E; 5 handled in detail in Sections 13.21 13.31 This 
completes the proof of Theorem [TJ 

Q.E.D. 



Note that, since the limit problem (|2~22|) . (|2~23| admits a unique solution for any square integrable initial 
data, there is no need to consider subsequences. 



4 Isotropic scaling: Proof of Theorem [2] 

We consider the situation where the Rossby and Mach numbers have the same scaling. Accordingly, 
we have m = 1, together with the uniform bounds derived in Paragraph 12.41 We recall that the static 
solution g e = g(\xh\) is now independent of e and we set 



Qe - Q 
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4.1 Preliminary results 

In accordance with the uniform bounds established in (| 2. 1711 2.15)) and (|2. 14112. lHj) , we may assume (up to 
taking subsequences) that 

q s -> g in L°°(Q,T;L~<(K)), r £ -> r weakly-(*) in L°°(0, T; L 2 (iT)) (4.1) 

for any compact if C f2, recalling that 7 > 3. Moreover, 

u e -> U weakly in i 2 (0, T; W 1 ' 3 ^; R 3 )). (4.2) 

The main point is to derive the equations satisfied by r and U. Similarly to Section [3.21 we rewrite the 
Navicr-Stokcs system in the following form: 

ed t r £ + div x (g £ u £ ) = 0, 

ed t (g £ u £ ) + ediv x (g £ u £ ® u E ) + V x (p'(g)r e ) + b x {g E u e ) ^ ^ 

= -£div x S(V x u e ) - -V, (p(g e ) - P {g) - P \g){g e - g)) + r e V x G. 

£ 

The bounds (|2. 14112. 15} allow to pass to the limit in the previous system, in the sense of distributions, to 
derive the following diagnostic equations: 

div x (£U) = 0, V x (p'{g) r) + bxgU = rV,G. 

Moreover, introducing 

R = P'(g)r, 

one can write the previous equations as 

div x (£U) = 0, V x i? + bxU = 0. 

From there, one has easily that 

R = R(x h ), U = [U h (x fc ),0] (4.4) 

together with 

div h (§U h ) = 0, V h R + Ui = 0. (4.5) 

Note that in addition to the constraint div^ (gUh) = 0, the last equation implies that div^U^ = 0, 
and V x g ■ Uh = 0, therefore R = R(\xh\) is radially symmetric. 

These constraints do not determine the time evolution of r and U^. Our aim is to show that, in accordance 
with the conclusion of Theorem [5J the function r satisfies equation (|2.26|) . 

The rest of the paper is devoted to the proof of (|2.26[) . To this end, we perform the analysis of the high 
frequency waves generated by the linear part of system (|4.3|) . Indeed, one can recast (|4.3|) in the general 
form 

ed t U e + L[U e ] = eN e + eT £ 
where U e = [r e , V e ], with V e = g £ u e , while 

L[U e ] = [div x V £ , gW x (P'(g) r £ ) + bxV E ] 
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stands for the linear part of the equation, and 

N E = [0, N e ], N e = -div x (g E u E ® u e ) . (4.6) 
Furthermore, T E includes the remaining terms: 

JT £ = (0,F £ ), F e = -diy x S(V x u e )-^V x (p(g e )-p{g)-p'(g)(g E -g)). (4.7) 

It can be checked that L defines a skew-symmetric operator with respect to the scalar product 

(U,W) = J (rr'P'(g) + V-V'f 1 ) dx; 

whence it generates time oscillations with frequency e , whose nonlinear interaction is not obvious. In 
fact, the the principal difficulty of this part of the paper is the fact that the spectral properties of the 
operator L are unclear, and thus, in contrast with the previous part, employing dispersive effects seems 
difficult if not impossible. We shall therefore rely on local methods and use as much as possible the 
structure of the equations, to control the asymptotic behavior of the convective term. 

4.2 Handling the convective term 

The goal of this section is to better understand the behavior of N e , and more precisely of 

N £ = div h Q (V E:h <g> V e>h )J , where V B>h = g e u e>h . (4.8) 

Later on, we shall compare this term with the average with respect to X3 of 

-N c ,/, = div^ (g E u E ® u e>h ) , 

which reads 

(-Ne,/,) = di V/l ({ggUe.fc (8 Ue,/,)) . (4.9) 

These are two main steps in deriving the limit equation (|2.26j) . We start by showing that 

(N B |VjJ-v) -> 0, e^O, 

for all -0 € ^^((O^T) x R 2 ) such that V/,g- V^V = 0. Here and hereafter, the symbol (-|-) denotes the 
standard duality pairing. We follow a general approach initiated by Lions and Masmoudi |17j . |18) in 
the framework of the incompressible limit, and later adapted in jllj to the case of rotating fluids. This 
approach, used also in the last step of the proof of Theorem[TJ is reminiscent of compensated compactness 
arguments. It relies on various cancellations, obtained by multiple use of the structure of (|4.3[) which 
roughly speaking reads 

ed t r £ + div x (g e u E ) = 0, 

ed t (g E u E ) + gW x (P'(g)r £ ) + b x (g E u E ) = 0(e) 
and will provide some compactness for appropriate combinations of r E , g E ,u E . 
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Similarly to Section I3~2l the complete treatment of N e involves spatial regularizations: 

u £ .5 = K5*u e , V £ ,a = Ks*(g e u e ), r £i s = n s *r £l 

where Kg = Kg(x) is a family of regularizing kernels. Note that the quantities V 6 g and r £t g are bounded 
in L^^T-L 1 ^)) for any compact K C Vt, and u M in i 2 (0, T; M /1 ' 2 (f2; i? 3 )) ', uniformly in <5,e. We 
claim the following result. 



Proposition 4.1 The fields u £ _g, *V £ ,$ and r £ j satisfy the following properties 

1 



V e , d = etlj + tl s , and curl, ( - V e , s ) = eT* s + T 2 4 



where 



\\' t l,s\\L 2 (0,T : W k ' 2 (K:R s )) + \\Tl,s\\L 2 (0,T:W k ' 2 (K:R s )) < c(5), 



||t 2 5 ||l2( T;H /i,2 (if< - ;J? 3)) + ||T 2 5 || L 2((o^T)xK;_R3) < C 

/or any compact if C and fc = 0, 1, . . . Moreover, 

sup \\r £ — r £i s\\L°°(o t-w s ' 2 (K)) ^ as 5 —> /or aZZ s < 0, 
sup ||u e — VL e ,s\\L 3 (o,T;W*< 3 (K;R 3 )) ^ as o" — > /or aZZ s < 1, 

Sup ||V e - V e,s\\L 2 (0,T;W>< 2 (K;R 3 )) -> as <5 ~> for all S < -1/2 

/or any compact K C f2. Finally the following approximate wave equation is satisfied: 

sd t r £t s + div x V £tS = 0, 

e9 t V £ , 5 + £V X (P'(£) r e ,,) + b x V B ,« = eF* 4 + F 2 ^, 
where the source terms are smooth and satisfy 



(4.10) 



(4.11) 



(4.12) 



™P\\Fl. S \\m0,T;WX.yK;R3)) < C{6) 



e>0 



SU P ||Fe,«|U=>((0 > T)xJr;fl3) + 

£>o y 

/or any compact K C f2, fc = 0, 1, . . . 



curl^-F^) 



as 5 



(4.13) 
(4.14) 

(4.15) 
(4.16) 



L 2 ((0,T)x_ff;iJ3) 



We postpone the proof of this proposition to Section 14.31 This regularization process allows to establish 
the following convergence result. 



Proposition 4.2 The nonlinear quantity 



satisfies 



N £i5 = div h - (V BiS ,h ® Ve,6,h) 
\Q 



(4.17) 



lim limsup 

s-yo y 



/or aZZ ^ G C* c °°((0,T) x i? 2 ) sucft ttai V h g ■ V^V = 0. 
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The rest of this part is devoted to the proof of Proposition 14.21 Actually we shall prove more precisely 
that 

N £ ,«5 = F eiS V h (j) + g £ .s^hQ + S e ,5 

for explicit functions 4> = <i>(t,%h), Fe,8 = Fe,s{o), 9e,s = 9e,s{t,x h ) and a remainder s £i5 = s e ^(t,x h ) 
satisfying 

lim f limsup | (s e s \<p) I ) = for all tp e C?°((0,T) x fl;R 2 ). (4.18) 

•5^0 \ £ ^ ' J 

This will of course imply the result. 

The notation o(l) will refer hereafter to any term s e _s satisfying (|4.18[) . Moreover we drop from now 
on the lowerscript 5, except if some ambiguity is liable to occur. 

Let us set 

N £ = Nl + N 2 £ = divfc Q (V e ,h) <g> (V £ , h )) + div h Q ({V e , h } <8> {V £jl }) 

where the notation (•), {•} has been introduced in p.26[) . and treat these parts separately. 

• Treatment of N*: 
We notice that 

N* = dw h ( (V 8)A ) ® 4 (V E . A )) = 4 (V e , fc ) di V/l ( (V e , h » + (V Blh ) • v h f i (V e>fc ) 

= i(V e , h >div fc ((V e>h » + lQVh(\l(Ve,h)\ 2 ) + (cml h l{V e , h ))(V e , h ) X . 
Q 2 V q J Q 

On the one hand, averaging (|4.13[) with respect to X3, and multiplying by (V £i /j), we get that 

\ (V e ,h) div h ((V e>fc )) = -e(ft (r e ))i (V e , fc ) = -eft (4 (r e ) (V^)) + 4 (r £ ) e (ft (V e>fc » , 

where the first term at the right-hand side is of order o(l). As regards the second term, averaging the 
horizontal components of (|4.14p with respect to £3 and multiplying by (r s ), we end up with 

\ (r e ) e (ft (V,.,)) = -V„ (P'(g) (r e » (r e ) - 4 (V e>fe ) x (r E ) + 4 < e F* fc + F^) (r e ) 

£> Q Q 



P^ V "Ql P '^)^)i 2 ) - |<V B) /.> X (r e ) + o(l), 



where the o(l) comes from the properties of the F* 's and the fact that r e = r £ g is uniformly bounded in 
Lf oc ((0,T) x il) with respect to e and (5, see (|4.1[) . Thus N* can be written in the form 

N* = ffcurlfci (V £)fe )) - 4 (r £ )) (V £ , h ) X + V ft + o(l). 
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It remains to handle the first term at the r.h.s. Therefore, we average (|4.14[) with respect to x$ , divide 
by g, take the curl of the horizontal components, and subtract average of (|4. 13[1 with respect to X3. We 
obtain, 

ed t ((cm\ h l (V e , h )) - 4 (r £ )) = 4 (V E ,h) ■ V h g + F E h (4.19) 



where 



F^curl^sF^ + F 2 ,). 



According to Proposition 14. II there is a function f(S) going to zero with S such that, uniformly in e 

1 



sup 

£>0 



cm\ h -F z h 

Q 



< /(<*)■ 



(4.20) 



L 2 ((0,T)x/f;i?2) 



Now let us notice that by definition of g one has 

^hQ P'{q) = 2x h , so V h g p'(g) = 2gx h . 
By Assumption (|2.2[) . this gives in particular that there is a constant C such that for all x^, 

\VhQ(x h )\>C\x h \. 
Now let us consider a smooth function \s defined by 



(4.21) 



Xs(xh) ■= X 



where x is a function of C^°(i? 2 ; [0, 1]) such that x( x h) = 1 if \%h\ < 1- Using (|4.10|) we get 



X s (curU- (V e /,)J - - (r, 



Vs./O = X* (<3 + ^3 ~ | (r.)j « fc + t% h y 
Now we can write, by Holder's inequality and the continuous embedding W 1,2 C L 6 



Xs ( < 3 + T £ 2 3 - - (r e ) ) « fc + t 2 ^)- 



< IIx5||l3(a-) 



L 2 (0,T:L 2 (K)) 



L 1 ((0,T)xif;_R 2 ) 



Recalling that r e is uniformly bounded in e,5 in L°°(0, T; X 2 oc (f2)), and using (|4.11|) and f|4.21|) we infer 
that 



X5 ( eTl 3 + T 2 3 - (r e ) ) (etl ih + t 2 h )~ 



< (ere 



L 1 ((0,T)xK;i? 2 ) 



(*) + c)|{s hj |V A e(ar fc )| < VM} 



< (ec(S)+c) B(0,Cy/f(S)) 

= 0(1). 
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This allows to conclude that 



Xs 



) x = o(l). 



Then we can write 



(1 - Xs) (( curi 4 (V e , h >) - | (r £ >) (V eih )" 



(1-Xtf) curl^(V E , ft ) )--(r 8 ) 



(1-Xtf) ( (curlA|(V eih ; 



Q 



(V 



V |Vfcg| 2 
(V Bth ) • V h £ 



(V 



iVfcgl 



n|2 



|V h e| 



for 







(curhj (V e , fc >) 


» 







(Vg) 1 ■ Vhg 



As regards the first term at the right-hand side of the last equality, we use (|4.f 91) to derive 



(1-Xs) (curl fc -(V e K))--(r e ) 



= — £<9t 



(1 - Xi) 

2|v fc el a 



e (curl,^ (V e , fc > ) - <r e ) ) ) + £ 2 ( < 3 + ^3 - ^ (r e > ) ^^4(1 - X.) 



with the notation introduced in (|4.10[) . We get, by (|4.f lj) . 



< 



L 2 (0,T;L 1 (if)) 



|V h £| 



|V h e| s 



L 2 ((0,T)xif) 



v/w 



which is o(l) thanks to (|4.20p . Combining the previous inequalities leads to 

= g 1 VhQ + F(g)V h </> + o(l), 

as expected. 

• Treatment of N 2 : 

This time, we consider 
1 



= div h ( - ({V £ ,4 <g> {V £ . h }) 



\ {V e , fc }div fc ({V B , h })\ + l -g /v fc 



- {Vs.h} 
Q 



curl h (^{V £ , h }) {V e , h }- 



23 



As — {V £ h} has zero vertical average, we can write 
Q 



curl^iVe} = ( 9a " B ' h ), n e , h := {V^ - ^M"* {V e>3 } , lo £ := cml h -_ {V £>h } , 

g \ UJ e J g g g 

where we have set 9^~ 1 a = I {a), see (|3.26[) . Applying curl^ to the momentum equation (|4.14[) yields 

ed t n £ . h = -„{V e , h ) + |fl3 1 (curl a ,i(eFi + Bt)) h |, (4.22) 
ed t u s = -div fc i {V £ , fc } + fcurl fc i(eP^+P?) fc \. (4.23) 



From there, we deduce that 



/curl ft (~{V E)/t }) {V £ ,,J^ = (w 8 {Ve./i^) = g(w B eftn£ h ) + o(l) 

= -g((sd t uj e )Sl^ h ) + (ed t (guj E Sl^ h )) + o(l) 

= e/n^div h |{v eih }\ + o(i) 

Here, we have used repeatedly equations (|4.22j) . (|4.23|) . The error terms generated by 

(ids 1 curls* (fFl + F 2 e 

are responsible for the o(l) term, as can be verified using as previously ()4.10[) - (|4.11|) . Therefore 
= {V Sth }div h ({V £ , h })\ + g(n^ h divA{V e . h }^ + F(g)W h i> + o(l) 
= \ <({V e , h } + gn^) div fc ({V e , fc })> + /gni h {V £ , h }-V h ±\ + F(g)V,^ + o(l). 



By straightforward manipulations we have 



i({v £ ,4 + ^)a 3 {^3} = %Qm,3}({v 6lh } + en^)) - ^ v »l 

We can replace the first term by using equation (|4,14[) 



{V e , 3 } I' - I {K, 3 } iVfc- 



\ ({{Ve,h} + ^) div, {V £ }) = -- {{ed t {r £ }) ({V e ,4 + §n$ h )) 
which leaves us with 

N e 2 = ~{{ed t {r £ }){{Y £ , h }+~Q^ h )) + (§fl$ h {V £ , h } ■ V h ~\ + F(g)V h iP + g 2 V h ~g + o(l) 
= ({r £ } ed t ({Ve,h} + m £ , h )) + (§nth { V ch} ■ + F (e)V>^ + g 2 V h ~g + o(l). 
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Now, thanks to pTTIj) and (|T22l , 



~({r s }sd t ({V £ , h } + Qn^ h )) = -^—W h \P'(g){r e }\ 2 + o(l), 

so that resumes to 

N2 = (^^[Ve.fc] • V h |\ + F(g)V h ip + 9 2 ^h~Q + o(l). 

Finally, we proceed as previously in the case of N*, with the first term at the right-hand side, this time 
omitting for simplicity the cut-off near V^g = 0. We write 

({v £)fc } • v^~) n£ fc = ({v £ ,4 • v fc g) f ^'1% + n fe'J a fcg v fc gl 



= g ( e fi,n e , fc • V h g) "^'.^ V^g + g 3 V h ~g + o(l) 
where we have used (|4.22l) in the passage from the second to the third line. This yields 

({v Slfc }-V/ig)n£ fc = 5 Ve + o(i) 

and finally we get an expression of the form: 

N 2 = F{g)V h ip + g\7 h g + (1) 
as expected. Proposition [472] is proved. 



Q.E.D. 



4.3 Regularization process 

In this section we shall prove Proposition UTT] We start by recalling that 

ed t r e + div x (V e ) =0, 

and 

ed t V £ + ~gV x (P'(g) r s ) + b x V e = eN £ + sF £1 

with notation (|4.6|) and (|4.7[) . The first step consists in establishing some bounds for N £ and F e . 
The energy bound clearly implies that N e is bounded in L°°(0, T; W _1,1 (if; R 3 ))- As for F e , one has 

clearly that div K S(V x u e ) is bounded in L 2 (0, T; W~ l ' 2 (K; R 3 )), and 4jV K (p(&) - p(g) - p'(g)(g e - g)) 

is bounded in L°°(0, T; W~ l,1 {K; R 3 )). Therefore in particular 

N e + S e is bounded in L 2 (0, T; W~ 5 ^ 2 ' 2 (K; R 3 )) (4.24) 

for any compact K C CI. 
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Now let us proceed to the regularization. First we notice that 

V e ,5 = sk s * (r e u £ ) + k s * (gu £ ) =■ etl s + t 2 eS 

and 

cur^ Q V £i a^ = ecurlr (^k s * (r e u e )^ + curl^ Qk 5 * {gu £ ^J =; eT^ s + T 2 S 
so thanks to the L 2 bound on r £ and the W 1 ' 2 bound on u £ , wc deduce easily that for all k, K 

llte^llL^O.TWk^ifift 3 )) + H T £,jL 2 (0,7W fc ' 2 (if;fl 3 )) < c {&)i 
Ute^lU^O.TiWi.^A';/? 3 )) + |Te,ii|L 2 ((0 J)xA';Ji 3 ) < C 

uniformly in e (and i5 for the second bound). This proves (|4.11[) . The uniform bounds derived previously 
also give directly the convergences (|4.12p . 

Now let us turn to the wave equations. By convolution we get (with obvious notation) 

ed t r £}S + div K V M = 0, (4.25) 

and 

ed t \ e , s + Q V, (P'(g) r e , s ) + b x V e>d = eF l £>s + F 2 S (4.26) 

with 

F* 4 = N £/5 + F £ . s 

and 

Fig = gV x (P'(g)r^) - {gV x (P'(g)r e )) * k s . 
Clearly (|4.24p implies (|4.15p . so let us turn to the statement (|4.16|) . 

In order to see l|4.16p . we use [THl Proposition 4.1] (which forces the restriction 7 > 3) on compactness of 
solutions to P~2"5]l . (g^Hl), namely, 

\\r s — r £t s\\Lp(lo,T}-L 2 (K)) — > as J — > for any compact K C £1 and any p > 1, 

together with Lemma 3.3 (2) of [19]. Note that, compared with the situation treated in [T9J Proposition 
4.1], the present system contains an extra term b x V £ ,a already known to be compact with respect to 
the space variable. 

The vanishing of F 2 s (uniformly in e) follows directly. To handle the convergence of curl^iF^ s , we then 
notice that 

curl^F^ = A{g)Vr £ , s - (A(g)Vr £ ) * n s 

for some smooth matrix function A. Still using Lemma 3.3 (2) in |19) . we obtain the vanishing of 
curLF^ s , which together with the one of F^ s completes the proof of Proposition 14. II 

Q.E.D. 
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4.4 Conclusion 



Thanks to Proposition 14.11 we can conclude the proof of Theorem [3J We keep the notation N e N e 
and N e , a of Section SjU sec ifCS)) . flH|) , (|3T7|) . 

Let V = € C£°((0,T) x r2 ) such that v ^ ' VftV = 0. We write 



((-N B , h > |V^V) - (N E , 5 |V^V 



p e u ej/l <g> u £ ,fc - -V c ,i,h (8) V E ,6,h ) : Vfc (g) V h ipdxdt 

JQ \ Q 



< 



I I ((Vs,h - Ve,6,h) ® u e,fc) : (Vfc ® V^V) d.Tdt 

Jo Jo 



/ T / (Ve,6,h ® " : (Vfc ® Vjty) dxdt 



JO 



V 



e,,5,fc 



Vs ' h VeAh 1 - (Vfc <g> V^dxdt 



4,5 + 4,5 + 4 3 <5 



We have: 



V e ,fc • ((V/, (8 V^)u e ,fc - k* * ((Vfc ® V^V)««,h)) 



0(<J) 



'0 Jo 

uniformly in e, using that 

|| (Vfc (8 V/J"V)Ue,fc - K S * ((Vfc <8> V h L V)Ue,fc)||L2((0,T)xO;i?2) < C5\\u e ,h\\L^(0,T;W^(K;R 2 )) 

for some compact K containing the support of tp. Then, noticing that 

V £ TV 



one obtains easily 



Finally, we remark that 



— = u e + e— u e 

Q Q 



lis < c(S)e. 



V e - V £ j = (us * (f?u e ) - pu e ) + e (r e u e - « 5 * (r £ u e )) 



to obtain 



4 3 5 < c<5 + c(<J)e. 



Putting these inequalities altogether yields 



limsup limsup (J £ 5 + / E s + I* s ) = 0. 
5->o e^O 
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Combining this with Proposition l4.il we deduce that 

lim (N e , h |Vjty) = 

for all ip = ip(t,Xh) £ C%°((0, T) x R 2 ) such that V hQ • V^?/> = 0, meaning the function ip is radially 
symmetric. 

We are now at the point of getting the equation satisfied by r, u, cf. (|2.24[) . The horizontal part of the 
momentum equation reads 

dt{g e u E h ) + div x (g 6 u 6 ® u eJl ) + e~ 2 V^p(fc) + e~ 1 g E u^ h = [div x S(S7 x u e )] h + % V h G. 

We recall that V^G = V^P(p) = P'(g)S/hQ- We integrate with respect to £3 the last equation, and 
apply curl/j. We obtain 

d t cml h (g e u eJl ) + curl/jdiv^ (g E u E . h ® u E) h) + e _1 divh (g E u eJl ) 
= curl h div h Sh,h{Vh (Ue,h)) + curU ( — V/ t g 

where 

2 

Sft.fcCVfcUft) = n(VhVLh + V^u ft - -div fc (u/,)I fcih ) with I hJl identity matrix in R 2 . 

Continuity equation yields div^ (g e u e .^) = ~sd t (r £ ); we employ this fact and (|4.27[) . where we use a 
radially symmetric test function ip £ C^°((0, T) x i? 2 ) to get 

(d t (cwl h (g s u Sth ) - (r e ) ) - curl/ l div fc S/ l) / l (V ;E (u £i / l })|V'J 

£ JO JB? 

The first (convective) term at the r.h.s. goes to zero, whereas the second one is identically zero by the 
properties of ip. All other quantities converge easily to yield 

(d t (cvulh (gV h ) - (r) J - n^ h cm\ h V h ,i^\ = 

for any radially symmetric ip. Using the properties (|2.24| - [2~23|) . we arrive at (|2.26[) . 

Finally, repeating the same procedure with ip £ C^°([0,T) x R 2 ), ip radially symmetric, we obtain 

- {{cur\ h (gu h ) - r)\d t ip) - /it (A/ l curl/ s U/ s |i/') = - (curl/ l (^u ,/ l ) - (r )\i>\ t=0 ) , (4.28) 

where (uo, ro) are the weak limits of the family of initial data uo. E , ro.e- This justifies the initial condition 
stated in (|2.27l) . Moreover, in view of hypothesis (|2.8[) . we have 

VP'(g){ro) G L 2 {R 2 ), y/§(uo, h ) £ L 2 (R 2 ;R 2 ). (4.29) 
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Under these circumstances, it is easy to show that (|2.24[) - (|2.25|) - (|4.28[) admits a unique solution. 
Indeed, taking ip = P'(g)r(0) as a test function in f|2 . 2T[) . we check that (|4.29[) implies 

J [g\W h (P'(g)r(0))f + P'(g)\r(0)\ 2 ) dx < 

whence uniqueness follows from standard energy arguments. See [9j section 5.2]. Thus, there is exactly 
one accumulation point for the sequence (r E , u £ ), and the whole sequence converges to it. 
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